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ABSTRACT 

We present continuum observations of the infrared dark cloud (IRDC) G48.660.22 (G48) obtained 
with Herschel, Spitzer, and APEX, in addition to several molecular line observations. The Herschel 
maps are used to derive temperature and column density maps of G48 using a model based on a 
modified blackbody. We find that G48 has a relatively simple structure and is relatively isolated; 
thus this IRDC provides an excellent target to study the collapse and fragmentation of a filamen- 
tary structure in the absence of complicating factors such as strong external feedback. The derived 
temperature structure of G48 is clearly non-isothermal from cloud to core scale. The column density 
peaks are spatially coincident with the lowest temperatures (~ 17.5 K) in G48. A total cloud mass of 
~ 39OM is derived from the column density maps. By comparing the luminosity-to-mass ratio of 13 
point sources detected in the Herschel /P ACS bands to evolutionary models, we find that two cores 
are likely to evolve into high-mass stars (M* > 8M ). The derived mean projected separation of 
point sources is smaller than in other IRDCs but in good agreement with theoretical predications for 
cylindrical collapse. We detect several molecular species such as CO, HCO+, HCN, HNC and N 2 H + . 
CO is depleted by a factor of ~ 3.5 compared to the expected interstellar abundance, from which we 
conclude that CO freezes out in the central region. Furthermore, the molecular clumps, associated 
with the sub-millimeter peaks in G48, appear to be gravitationally unbound or just pressure confined. 
The analysis of critical line masses in G48 show that the entire filament is collapsing, overcoming any 
internal support. 

Subject headings: Stars: formation, Stars: massive, Stars: kinematics and dynamics, ISM: clouds, 
ISM: molecules, ISM: individual: G048.66-0.29 



1. INTRODUCTION 

Our knowledge of the physical conditions during the 
very early phases of high-mass star format ion (> 8M^) 
has b een relatively limited in the past (e.g., |Beuther et ah] 
2007), due to significant observational challenges. Con- 



sequently, the qualitative and quantitative differences in 
star-formation processes o ver different mass regime s are 
still strongly debated (cf. |Zinnecker fc Yorke||2007[ ). In 
the endeavor to enhance our understanding ol star for- 
mation beyond the well-studied low mass stellar nurs- 
eries, infrared dark clouds (IRDCs) have become more 
and more important targets. Initially IRDCs were dis- 



covered by the Infrared Space Observatory (ISO, Perault 
et al.|1996|) and th e Midcourse Space Experiment (MSX, 



Egan et al. 1998). They appear as filamentary opaque 
structures against the diffuse galactic mid-infrared back- 
ground. As IRDCs are known to be cold (< 20 K), dense, 
and exhibiting high column densities (> 10 23 cm -2 ), they 
are suitable environments to form young stellar objects 
that c an evolve into stars of intermediate to high mass 
(e.g., |Carey et al. 2000 Rathborne et al.| 2006; Ragan 
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x Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA, 



z Based on observations carried out with the 1RAM 30m Tele- 



scope. IRAM is supported by INSU/CNRS (France), MPG (Ger- 



many) and 1(J1N (Spam). 



|et al ||2006l [2009] |Vasyunina elHuT||2W)9| . IRDCs were 
widely discussed in the context ol massive star forma- 
tion, but a certain fraction of IRDCs only host sub- solar 
to intermediate-mass star formation (|Beuther et al. 2007 
iKauffmann fc Pil lai 2010; |Peretto fc^uller||2010| ).^tm 
|Rathborne et al.| ( |2006|) argue that most stars in our 
galaxy could lorm m IKDC-like molecular structures. 

Most of the current theories about the formation of 
IRDCs treat the collapsing gaseous medium as isother- 
mal. With the Herschel observatory, we have the spatial 
resolving power, spectral coverage, sensitivity and dy- 
namic range to study the temperature and column den- 
sity structure of IRDCs in more detail. Already one of 
the early Herschel studies using spatially resolved SED 
fitting toward an I RDC indicated a non -isothermal tem- 
perature structur e (iPeretto et al.||2010| ). 

[Wilcock et al.| pOT^ T compared Herschel observations 
of IRDCs without embedded IR sources to radiative 
transfer models of strictly externally heated structures. 
They find resulting temperature profiles that go down 
to 11.5 K (sample average) in the center. Two further 
studies analyzed the temperature struc ture using SED- 
fitting and a modified blackbody model. Battersby et al. 
([2011 1) present large scale temperature maps of IRDCs. 
Beuther et al.J (|2011|) studied the IRDC18454 near the 



W43 mini-starburst and present dust temperature maps. 
Both studies found temperatures higher by more than 
5K compared to the IRDCs fromlWilcock et all (120121). 
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Fig. 1. — The upper panel presents a Herschel map of the galactic neigh borhood of G48. The do minant feature to the left of G48 is 
the Giant Molecular Cloud complex of W51. Two HII regions described bylKuchar & Banial \ 19941) are located above G48. Below G48 
the infrared bubble N100 appears. The map is composed of the PACS band at YU/im (blue), loU/im (green) and SPIRE at 250 jum (red) 
from the Hi-GAL field 48_0. The bottom right panel is a zoom into the whole IRDC. It uses the same Herschel bands, utilizing the EPOS 
observations. The bottom left panel shows the central part of G48 as an IRAC 3-color image (3.5, 4.8, 8.0/xm). The green band shows 
elongated lobes attributed to shock-excited H2, s ome times refereed to as "green and fuzzy" features "or extended green objects" (EGOs). 
The detected Herschel point-sources from Section [3^5] are indicated by circles in the right panels. All intensities are plotted with logarithmic 
scaling in all maps. 
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Further investigations have to clarify if the majority of 
IRDCs show similar minimum temperatures or if clumps 
with higher central temperatures exist. This classifica- 
tion is essential for the understanding of the formation 
process of massive stars. The comparison between dust 
and gas temperatures can shed light on the efficiency of 
the gas-dust coupling at densities typical for IRDCs. Fur- 
thermore, a well constrained dust temperature structure 
can provide critical input for chemical models of IRDCs, 
and in particular for the dust grain surface reactions. 

The goal of this paper is to accurately characterize the 
physical conditions of (high-mass) star formation in a rel- 
atively isolated IRDC, G048.66-0.29 (henceforth called 
"G48"). This case study is based on observations taken 
as part of the Earliest Phases of Star Form ation (EPoS) 
Herschel guaranteed time key program ( Launhardt et| 



al. 



20121 iRagan et al. 20121). We present sub-par sec 



scale dust temperature and column density maps and 
place these properties in context with the kinematic en- 
vironment, probed by various molecular gas tracers (CO, 
HCO+, N 2 H + , HCN, HNC). Herschel also enables us 
to characterize the embedded protostellar core popula- 
tion and, through the use of evolutionary models, try to 
understand the likely outcome of star formation in this 
cloud. 

The kin ematic distance to G48 was de termined to be 
- 2.7kpc flOrmel et al.|2005[ ) 2.5 kpc ( [Simon et al.|2006b ) 
or 2.6 kpc (this work) . G48 appears to be projected near 
the GMC W51 ( Figure [11 left), w hich is in fact twice as 



distant as G48 (Sato et al.||2010 , see also Section 3.2). 
G48 appears isolated; and therefore we can assume it is 
not strongly influenced by local heating sources. Fur- 
thermore G48 has a quite simple monofilament ary struc- 
ture. By comparing the observational results from this 
isolated IRDC to more complex IRDCs (c.f. Rag an et| 
aI1[20T2t , we can constrain the dominant effects at work 
Tor star formation in IRDCs. Simon et al.| ( |2006b ) es- 
timated from 13 CO observations an average density for 
H2 associated with G48 of ~ 3 • 10 3 cm -3 , a peak column 
density of 8. 8 • 10 21 cm" 2 and a LTE mass of - 580 M G 



Ormel et al. (2005) identified three different clumps with 
masses from 91 to 130 M using SCUBA and MSX obser- 
vations. By modeling the internal velocity structure of 
these clumps using HCO+ (3 - 2) and (4 - 3) single-dish 
data they found that the line- widths decrease outside the 
cloud, indicating a decrease in turbulent motions out- 
wards of the clumps. This suggests that the star forma- 
tion activity inside these clumps down to the core scale 
is driven by turbulence. 



van der 
the 



Wiel 



Shipman 
G48, 



^ (|2008j) identified 13 YSOs 
based on Spitzer mid-infrared 



m tne vicinity ol 
maps. The majority of these objects belo ng to "Stage 
I" , referri ng to the classification scheme of Robitaille et 



al. (2006), and have low to intermediate masses. From 



a com parison of [CI] and 13 CO data, |Ossenkopf et al. 



(2011) inferred the presence of a young photodissocia- 
tion region (PDR) ionized by the YSOs inside of G48. 

Throughout the literature the limits of sizes and 
masses for clumps and cores vary. In t his wo rk the spa- 
tial definition from |Bergin fc Tafalla] ( |2007[ ) is applied: 
cores range from 0.03 — 0.2 pc, while clumps extend over 
0.3 -3pc. 

Given the isolated nature and simple structure of G48, 



we use this source to study the collapse and fragmenta- 
tion of a filamentary IRDC. Through the characteriza- 
tion of the point sources the evolutionary stage can be 
assessed a final stellar mass can be estimated from differ- 
ent models. Furthermore, the Herschel maps are utilized 
to study the temperature and column density structure 
of G48. Molecular line observations provide information 
about the morphology of the clumps and clues to the 
pressure support therein. We assess the level of freeze- 
out by studying the CO-depletion in the central part of 
the IRDC. 

The next section describes the observations, the data 
reduction and some basic data processing steps for Her- 
schel, Spitzer, sub-millimeter continuum and molecular 
line observations. Section [3] details th e sc ientific analy- 
sis and the direct results. In Section 13.11 the structure 
and environment of G48 is described. The measured dis- 
tances for G48 and other sources in the near projected 
distance are discussed in Section |3.1| The temperature 
mapping, utilizing Herschel, SCUBA and ATLASGAL 
data, is delineated in Section [33] The subsequent Sec- 



tion 3.4 compares the column densities with other ob- 
serv ation s and estimates mass for the whole IRDC. Sec- 
tion |3.5| addresses the fitting and classification of point 
sou rces from the MIPS 24 /im and PACS bands. Section 
3.6 contains the analysis of the molecular line data. The 
concluding Section [4] discusses the results in respect to 
previous studies of IRDC and theoretical considerations. 
The last Section [5] summarizes the main results. 

2. OBSERVATIONS, DATA REDUCTION AND 
PROCESSING 

2.1. Herschel 

In this paper we present results from the Herschel 
GTO key program - "Early Phases of Star Formation" 
(EPoS). The high-mass portion of EPoS is compiled from 
different surveys. The sample includes a sub-sample of 
"classical" IRDCs, a selection of IRDCs in the near vicin- 
ity of high-mass proto-stellar objects (HMPOs) and cold, 
massive, compact sources detected by the ISO Serendip- 
ity Survey (ISOSS). For a detailed discussion of this sam- 
ple and the underlining selection strategy we refer the 
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Fig. 2. — The LL-spectrum from Source 2 from IRS/ Spitzer is 
shown. In the left panel the source spectrum is plotted in black 
and the extracted background spectra are shown in blue and green. 
The background spectrum extracted in the first nod position (BG 
Nod 1) pointed towards the GMC, while the second background 
spectrum is taken at the western edge of the IRDC (BG Nod 1). 
The atomic [SIII] line is present as a bright extended emission over 
the whole slit. The CO2 ice absorption band is underlined in light 
grey. The right plot shows the baseline subtracted CO2 profiles. 
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Fig. 3. — The intensity plots for G48 are shown in the different Spitzer, PACS and SPIRE channels. Each plot is shown on a linear scale 
to highlight absorption or emission at each w avelength band; the top color bar indicates the relative flux scale. The positions of the 24 /im 
point sources as characterized by van der Wiel & Shipman (2008) are indicated by red crosses in the MIPS 24 and PACS 70 /im maps. 
The MIPS 24 /im sources from van der Wiel &; Shipman (20U8) are labeled "Sx" while two additional 24 /im sources associated to G48 are 
labeled "El" and "E2". The detected sources m the JrA(JS bands are marked by the white circles. The radius of the circles imply the 
fitting radius between the different bands. The SPIRE maps are overlayed by SCUBA 450 /im, 850 fim and ATLASGAL 850 /im contours, 
respectively. The corresponding beam sizes fo r Syi tzer, PACS and SPIRE are indicated with circles in each panel. An isothermal contour 
of 19 K from the temperature fitting (Section |3.3| is overplotted for IRAC 8 /im and PACS 100 /im. The beam size of the observations 
plotted in contours are plotted as open ellipses. The white bar (right side) gives a projected scale of 0.5 pc on the map. 
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Fig. 4. — The integrated line intensities of the indicated species are shown as contours overlaying the 250 fim SPIRE map. The contours 
in the top and middle panels refer to the MOPRA observations, while the bottom panels present the CO data. The beam sizes are indicated 
in the last panel. The maps are integrated between 32 and 36km/s. The contours for the different line transitions are: 
HCO+: from 0.3 to LSKkms -1 in steps of 0.3Kkms _1 ; HCN: from 0.2 to O^Kkms -1 in steps of O.lKkms -1 ; 
HNC: from 0.2 to 0.8Kkms _1 in steps of O.lKkms -1 ; N 2 H+: from 0.5 to 2.0Kkms~ 1 in steps of 0.25Kkms _1 ; 
13 CO(1-0): from 0.8 to 2.9Kkms _1 in steps of 0.7Kkms _1 ; 13 CO(2-l): from 2.0 to Q-SKkms -1 in steps of 1.5Kkms _1 ; 
C 18 0(2-l): from 0.7 to 2.7Kkms _1 in steps of 0.5Kkms _1 . 
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reader to the overview paper by Ragan et al. (2012). 

The EPoS observation were supplement with data from 
the the "Herschel infrared Galactic Plane Survey" (Hi- 
GAL) project, done with the Photometer Array C am- 
era and Spectrometer (PACS, |Poglitsch et al.|2 010) and 
the Spectral and P hotometric Imaging Receiver (SPIRE, 
Griffin et aL|201 |) on board o f the Herschel Space Ob- 
servatory ( |Pilbratt et al.||2010[ ). 



2.1.1. PACS 

The EPoS PACS scan maps at 70 /im, 100 jam and 
160 /im were obtained with the medium scan speed of 
20" /s and the nominal scan leg length of 7'. The data 
were processed to level- 1, i.e. c alibrated and converted 
to physical units, using HIPE ( |Ott||2010| ). We applied 
a second level deglitching to flag bad data values by a- 
clipping flux values mapped on to each pixel. The 'time 
ordered' option was used, and a 25 a-threshold was ap- 
plied. The level-2 products , i.e. fully red uced fits maps , 
were obtained using SCANAMORPHOS ( |Roussel|[2012 ) 
version 5.0 instead of applying the nominal HIPE high- 
pass filtering and phot Project reduction step, due to the 
improved perfor mance of SCANAMORPHOS (see 



gan 



et al. 



Ra- 



2012). The maps in each band were reduced 



with a pixel-size of l'/O/px 

2.1.2. SPIRE 

The EPoS SPIRE scan maps at 250 /im, 350 /im and 
500 /im were obtained with a scan speed of 30" /s. They 
were reduced up to level- 1 using HIPE version 5.0with 
the calibration pool 5.1. The final level-2 reduction 
was done with SCANAMORPHOS version 9.0 using the 
'galactic' option. The pixel sizes are 6.0" /px at 250 /im, 
10.0" /px at 350 /im and 14.0" /px at 500 /im. 

2.1.3. Hi-GAL 

To include a larger coverage we obtained the raw data 
products from the Hi-GAL Open Time Key Project. 
These observations were carried out in the PACS/SPIRE 
parallel mode with a scan speed of 60"/s, resulting in 
large scale maps. The 70 jam and 160 jam PACS maps 
have been reduced using HIPE 8.0, applying the same 
options as for the EPoS PACS data. Level-2 process- 
ing has been performed with SCANAMORPHOS ver- 
sion 15. The pixel size is 3'/2/px for both maps. The 
250 /im, 350 jam and 500 jam SPIRE maps were reduced 
up to level- 1 using HIPE version 8.0 with the calibration 
tree 8.1. SCANAMORPHOS version 15 was used to pro- 
duce the level-2 maps. The pixel size is the same as for 
the EPoS SPIRE maps. The final maps are presented in 
the left panel of Figure [I] The Hi-GA L ma ps are used 
for the temperature mapping in Section pO 



2.1.4. Photometry, point-source detection 

Due to the diffuse and heterogeneous background emis- 
sion varying on a wide range of spatial scales, it is more 
robust to perform PSF photometry than aperture pho- 
tometry. Therefore the star finder PSF photometry 
software ( jDiolaiti et al.|[2000[ ) was used for the point 
source detection and flux extraction for the PACS bands. 
The latest empirical PSF of the Vesta asteroid was used 
(OD 160), obtained from the Herschel Science Archive 
(HSA), since an empirical PSF derived from sources in 



G48 is noise dominated and therefore less acceptable. For 
consistency the PSF maps were reprocessed with SCAN- 
MORPHOS and rotated to the orientation of the maps. 
The noise in the PACS maps is not white noise, since it 
suffers from pixel-to-pixel correlations induced, e.g., by 
non-thermal low- frequency noise. To calculate the aver- 
age noise in the map a Gaussian distributed noise was 
assumed, but corrected by a factor for the pixel-to-pixel 
correlation. 

The major source of confusion in the source detection 
is the complicated background. In particular the bright 
background towards the GMC returns numerous false de- 
tections. We only include regions of the PACS maps 
that cor resp ond to regions in the ATLASGAL map (see 
Section 2.3) with fluxes above a 3 a level. To further 
improve the detection efficiency an unsharp-mask-filter 
was applied in each PACS-band by subtracting a me- 
dian smoothed map from the respective science- frame. 
With this procedure we avoid false detections, and miss 
three obvious sources for the automated detection run. 
These sources were then included in the catalog identi- 
fied by eye as clear detections. With the compiled fixed 
source positions the final photometry was performed on 
the non-filtered science frames. To be considered a de- 
tection a source must meet the following conditions: it 
must have a source flux above the 5<r noise level and be 
detected in all PACS bands within the FWHM of the 
160 jam PACS PSF (~ ll'/2). The results of the point 
source extraction are reported in Section |3l?l The given 
fluxes are already color-corrected for the source tempera- 
ture^] Further details of the source extrac tion are given 
in th e EPOS high-mass overview paper (Ragan et al. 
20121. 



2.2. Spitzer 

We obtained the Spitzer IRAC and MIPS maps from 
the Galactic Legacy Infrared Mid-Plane Survey Extraor- 
dinaire (GLIMPSE, |Benjamin et al.|2003| a nd the MIPS 
Inner Galactic Plane Survey (MIPSGAL, |Carey et al.| 



2009] ). The 24 /am MIPSGAL maps were pro cessed witrT 
HieMIPS Data Analysis Tool (DAT v3.10 



Gordon et 



aT||2005[ ) starting from the raw data products. 

We used the 8 jam IRAC map to create an extinction 
map. PSF photometry was performed on this map, using 
the starfinder package. An averaged, empirical PSF 
was used to remove all the point sources in this field. The 
residual map was used to compute the optical depths. 
One of the key problem in calculating an optical depth in 
IR data, is to determine the foreground and background 
intensities. Sev eral different methods are described in the 
literature (e.ft |Simon et al.||2006al [Butler fc Tan||2009| 

However, to 
cloud we used 



Ragan et al. 

estimate the optical depth 



2009 



IButler fc Il5||2012| ). 

I depth in one single 



the method described in Vasyunina et al. (2009). 



In 



(.. 



Virdc _ h 



(1) 



where /irdc is the flux derived from the IRDC directly, 
Io is measured of "off-cloud" in a nearby region meant 
to represent the characteristic unattenuated MIR emis- 

3 See: http://herschel.esac.esa.int/twiki/pub/Public/ 
PacsCalibrationWeb/cc_report_vl.pdt 
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sion. Io (and /irdc) contains contributions from both 
the foreground and background (i.e. Iq = I\> g + if g ), and 
only Jbg is needed to computed the optical depth. For 
Jf g , we assume the zodiacal light emission provided in the 
header of the IRAC data product. This is a lower limit 
to the true foreground. Nevertheless, the uncertainties in 
the column densities are dominated by what is assumed 
about the extinction law (factor 2 to 3), whereas the if g 
contributes about an error less than 50%. 

T he p oint-source fluxes from MIPS 24 /xm used in Sec- 
tion [33] were also obtained with star finder. For MIPS 
PSF photometry we utilized a model PSF computed by 
the STINYTI# ] program, based on TINYTIM for the HST 
flKrist |[2QQ2| ). 

The spectrum, shown in Figure [2j was obtained from 
the Infrared Spectrograph (IRS) on-board the Spitzer 
Space Telescope. The observation (AOT: 12093696) was 
taken from the Spitzer Heritage Archive (SHA) as part 
of the "Infrared Dark Clouds in the Inner Galaxy" pro- 
gram. We used a pipeline partly based on the SM ART 



package and the FEPS spectral extraction tools flBouw - 



man et al. 2006 



Fhe background subtracted on- source spectrum and the 
background spectra from the off-source nod-positions is 
shown in Figure [2 Only the long wavelength-order (LL) 
is reduced, since the slit for the short wavelength did not 
cover any point source in the IRDC. The LL slit did only 
cover about 60% of the SMART-source PSF, which leads 
to a flux leakage in the continuum flux. But line shapes 
and line fluxes should be unchanged for the source spec- 
trum, except for a lower S/N-ratio. A strong [SIII] emis- 
sion appears over the whole IRS-slit, leading to an arti- 
fact in the source spectrum due to an over-subtraction 
of the background. The source spectrum was corrected 
for this artifact using a polynomia l background fit (see 
Section 3.2.3 in Pitann et al.||20lT ). 



2.3. Sub- Millimeter imaging 
We utilized sub-millimeter data from the SCUBA 



Legacy Catalogue (Di Francesco et al. 2008). The 450 



and 850 /im observations were part of the fundamental 
catalog with a effective beamsize of 17" 3 and 22"9, re- 
spectively. Given the better sensitivity of the survey, the 
870 /im observations from the "APEX Teles cope Large 
Ar ea Sur vey of the Galaxy" (ATLASGAL, |Schuller et 



ah 2009) were used in the temperature fitting ( Se ction 
X5) and to determine column densities (Section 3.4). 
The corresponding beam size for the ATLASGAL data 
is around 19'/2. The RMS noise level for ATLASGAL is 
~ 50 mJy. 

2.4. Millimeter line observations 

We obtained molecular line data using the 22 meter 
MOPRA radio telescope and the MOPRA spectrome- 
ter (MOPS), operated by the Australian National Facil- 
ity (ATNF). The observations were performed with the 
3 mm band receiver (90 GHz) in the "On-the-fly Map- 
ping" mode (OTF) with a beam size of 36" at 86 GHz 
and 33" at 115 GHz. Using MOPS in the "Zoom" setup 
allowed us to observe 13 lines simultaneously ( 13 CS, C2H, 

4 See: http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/tools/contributed/general/stmytim/ 



CH 3 CCH, CH 3 CN, H 13 CN, HCO+, H 13 CO+, HCCCN, 
HCN, HNC, HNCO, N 2 H+, SiO). 

Most of these lines a re good tracers for dense gas (see 
|Vasyunina et al.|[2011[ for a detailed description of each 
molecular species). However, the RMS noise for all MO- 
PRA maps is only slightly below 0.1K, that allows us to 
detect line emission above 3 a only in the HCN, HNC, 
HCO+ and N 2 H+ maps. 

The spectrometer configuration results in a velocity 
resolution of ~ 0.11 kms -1 . The main beam effic iency 



varies betw een 49% at 86 GHz and 44% at 100 GHz ( |Ladd 
|et al.|[2005| ). The typical system temperatures were be- 
tween 210 and 290 K. The 13 CO (1-0) data prod ucts were 



obtained fr om the Galactic Ring Survey (GRS, [Jackson 
et al.|2006 ) performed by the Five College Radio Astron- 



Swain et al. 2008 Pitann et al. pout), et al 



omy Observatory (FCRAO). With a angular resolution 
of 46" the beam of FCRAO is somewhat larger than that 
of MOPRA. The RMS noise level in the GRS map is be- 
low 0.9 K. 

We obtained 13 CO (2-1) and C 18 (2-1) m aps using 
the 230 GHz heterodyne receiver array (HERA, |Schuster| 



2004[ ) at the IRAM 30 m telescope. The average 
FWHMot the HERA beam is ll'/3, which is much better 
compared to MOPRA or FCRAO. The RMS noise level 
is below 0.5 and 0.3 K in the 13 CO and C 18 HERA 
maps, respectively. 

To acquire methylacetylene (CH 3 C=CH) data obser- 
vations were performed with the IRAM 30 m telescope 
in position switching mode. For these observations the 
EMIR receiver with Fast Fourier Transform Spectrome- 
ter (FTS) as a backend was used. Total time On source 
and Off source was 33 min. The Off source position was 
chosen 600" away from the On source position. The On 
source position was pointed at 19 h 21 m 49 s .9 +13°49 / 34 / /7 
(FK5, J2000). The pointing was checked every hour, 
giving a pointing accuracy better than 3". The typi- 
cal system temperature during observations was about 
200 K. At 153 GHz the telescope beam size is 16" and 
the beam efficiency is 0.73. The data were analyzed us- 
ing GILDAS/CLASS software. 

A list of physical parameters for each observed line is 
given m Table [1] 

3. RESULTS 

3.1. Morphology - From large to small scales 

The left panel of Figure [I] shows the galactic neighbor- 
hood of G48. The image is composed of the 70, 160 and 
250 jam Hi-GAL maps in blue, green and red, respec- 
tively. A global brightness gradient towards the outer 
parts of the galactic plane is present over the whole im- 
age. In the vicinity of G48 the GMC W51 is observed. 
The cloud complex W51A and W51B have an angular 
distance of ~ 15' and ~ 50' to G48, respectively. The 
EPoS maps (Figure fl] top right panel) , in particular the 
PACS observations, show a brightness gradient towards 



W51 and the infrared bubble N100 ( Churchwell et al. 
2006). W51 and G48 are not physically connected since 



recent observations place W51 twice as far away as G48 
(see next Section) 

In the Spitzer IRAC 8 /im and the MIPS 24 /im chan- 
nel the IRDC appears as the typical dark filamentary 
absorption feature against the diffuse galactic emission 
background (see Figure[3|. G48 is extended over ~ 7.' 5 in 
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TABLE l 

Molecular lines 



Molecule Trans. u r AV t h errn A u i E u /kB n cr it beam size Instrument Obs. Date 

GHz kms" 1 10 -5 s -1 K cm -3 arcsec 



HCO+ 


1 - 





89.188523 


0.166 


3.0 


4.28 


2.1 • 


10 5 


38" 


MOPRA 


2010 Jul 7-9 


HCN 


1 - 





88.631602 


0.172 


2.4 


4.25 


1.5 • 


10 6 


38" 


MOPRA 


2010 Jul 7-9 


HNC 


1 — 





90.663568 


0.172 


2.690 


4.35 


3.2 • 


10 5 


00// 
38 


MOPRA 


2010 Jul 7-9 


N2H+ 


1 — 





93.173700 


0.166 


3.6 


4.47 


1.8 • 


10 5 


00// 
38 


MOPRA 


2010 Jul 7-9 


CU 


1 — 





110.Z01354 


0.166 


r\ r\r\r' o 

0.0063 


5.29 


1.9 • 


10 3 


A All 

44 


b ORAO 


Archive data 1 




2 — 


1 


zz0.oyoo84 


0.166 


0.0b04 


15.87 


2.0 • 


10 4 


ll" 25 


TT) A TV /r /TJTT'O A 

1RAM / Hl^ltA 


20 12 Jan 26 


c 18 o 


2 - 


1 


219.560319 


0.163 


0.0601 


15.81 


1.9 • 


10 4 


ll'/25 


IRAM/HERA 


2012 Jan 23 


CHgc^H 


9 - 


80 


153.8172 


0.142 


8.75 








16" 


IRAM /EMIR-FTS 


2011 Jun 10/12 


CHgc^H 


9i - 


81 


153.8142 


0.142 


8.44 








16" 


IRAM /EMIR-FTS 


2011 Jun 10/12 


CHgc^H 


9 2 - 


82 


153.8054 


0.142 


7.55 








16" 


IRAM /EMIR-FTS 


2011 Jun 10/12 




9 3 - 


83 


153.7907 


0.142 


12.52 








16" 


IRAM /EMIR-FTS 


2011 Jun 10/12 


CH3C2H 


5i - 


4i 


85.455665 


0.142 


0.1948 


19.53 






38" 


MOPRA 


2010 Jul 7-9 


H 13 CN 


1 - 





86.339860 


0.169 


2.2256 


4.14 


2.6 • 


10 6 


38" 


MOPRA 


2010 Jul 7-9 


H 13 CO+ 


1 - 





86.754288 


0.164 


3.8534 


4.16 


1.9 • 


10 5 


38" 


MOPRA 


2010 Jul 7-9 


SiO 


2 - 


1 


86.846960 


0.135 


2.927 


6.25 


1.9 • 


10 6 


38" 


MOPRA 


2010 Jul 7-9 


C 2 H 


2 - 


1 


87.316925 


0.179 


0.1528 


4.19 


2 • 10 5 ( a ' b ) 


38" 


MOPRA 


2010 Jul 7-9 


HNCO 


4 ,4 - 


3o,3 


87.925238 


0.137 


0.9024 


10.55 


4.5 • 


10 6 


38" 


MOPRA 


2010 Jul 7-9 


HC 3 N 


10 - 


- 9 


90.978865 


0.125 


5.812 


24.01 


7.5 • 


10 4 


38" 


MOPRA 


2010 Jul 7-9 


CH 3 CN 


5i - 


4i 


91.985316 


0.140 


6.0800 


20.39 


4- 1C 


5 (b) 


38" 


MOPRA 


2010 Jul 7-9 


13 cs 


2 - 


1 


92.494308 


0.134 


1.4124 


6.66 


5.0 • 


10 5 


38" 


MOPRA 


2010 Jul 7-9 



Note. — The columns list, the molecular species, the transition, the rest frame frequency, the thermal line width at 17. 3K, the Einstein 
coefficient, the upperstate energy, the critical densities and the angular resolution for our observations, the used instrument, from left to right, 
respectively. The upper part of the table presents the detected molecular lines, while the bottom part lists the MOPRA lines lacking a 3 cr-detection. 

1 from the Galactic Ri ng Survey with FCRAO 
Lo et al.N2009> 



from 
from 



Sanhueza et al. 



( 2012| 



the north-south direction. The dark cloud can be roughly 
divided into three parts. The central region of the IRDC 
appears as a elongated structure of ~ 2: [ 3 in east- western 
orientation. This region has the brightest sub-millimeter 
emission. The western part of the central region shows 
the highest extinction in the Spitzer maps. Only this 
part of the IRDC appears in clear absorption in all three 
PACS channels. Furthermore this part is clearly associ- 
ated with the most massive sub- millimeter clump, called 
"extinction peak" (EP, with a mass of ~ 130 M from 
SCUBA observation) in |Ormel et al.| (|2005). The eastern 
part of the central IRDC only shows clear extinction in 
the Spitzer bands. In the Herschel PACS maps this re- 
gion is dominated by compact background emission from 
emb edde d point sources, in particular Source 1 (see Sec- 
tion 3.5). The central region, as the rest of G48, appears 
as a amuse emission structure in the SPIRE maps, but 
two compact sources can be identified in the eastern cen- 
tral part. The peak positions of these emission features 
detected at 250/im and 350/im are in good agreement 
with the peak positions of two sub- millimeter clump s (PI 



with 91 Mq and P2 with 100 M ) identified by Ormel 
et al. (2005). Due to the lower resolution of the 500 /im 



SPIRE channel the two positions are not resolved. 

The structures found south- and northwards of the cen- 
tral part of the IRDC appear more diffuse in the Spitzer 
bands and are barely recognizable in the PACS chan- 
nels. In contrast, in the SPIRE bands this diffuse absorp- 
tion features appear as an elongated emission structure. 
While this emission is diffuse in the northern part, com- 
pact emission is detected in good agreement with the 
SCUBA and ATLASGAL data in the southern part of 
the IRD C. The associated sub-mill i meter peak was called 
"P3" in |van der Wiel fc Shipman] ( |2008[ ) (see Figure [3). 



3.2. Distance 

The LSR velocities found for the different clumps and 
molecular species do not differ by more than 3 ay . This 
results in a weighted me an v elocity of (Vlsr) — 33.91 ± 
0.05 kms -1 (see Section 3.6). The resulting kinematic 
distance for G48 is 2. 6 ±0.6 kpc. This value is used for all 
calculatio ns in this paper. W e derive kinematic distances 
using the Reid et al. (20091) algorithm, which takes into 



account trigonometric parallax measurements of masers 
in high-mass star forming regions in the Galaxy, applied 
to a disk Galaxy model. We use their standard derived 
rotational parameters of 0o = 254 km s _1 and Ro = 
8.4 kpc. 

The Hi- GAL three color image (Figure pi) shows G48 at 
a close projected distance to the GMC W51 and several 
objects, as e.g., HII regions, infrared bubble and other 
bright point-sources. Water m aser observation e stimates 
a distance of ~5.4kpc for W51 ( Sato et al.|2010 ). Known 
distances of infrared bubbles (Sarah Kendrew 2012, pri- 
vate communication) and sources from t he RMS survey 
catalog (e.g., massive YSO, HII regions |Urquhart et al.| 
2008 ) in the vicinity of G48 have a distance ol more than 
5 kpc. Moreover, no CO is detected in the global GRS 
maps outside of G48 at its typical LSR velocities. 

We therefore conclude that G48 is an isolated region 
making it an unique IRDC. 

3.3. Herschel temperature map 

To create temperature and column density maps of 
G48 we used Herschel maps and ground based sub- 
millimeter data. In brief, the basic steps we have taken 
to fit the temperature and column density maps are as 
follows: 
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Fig. 5. — The top two panels show H 2 column density maps obtained from IRAC 8jum extinction (left) and ATLASGAL 870 fim (right) 
observations. In these top panels the column density thresholds for the mass determination are shown as contours: within the contour of 
1.34 • 10 22 cm -2 a mass of 450 M is enclosed for IRAC 8 fim while the 0.85 • 10 22 cm -2 ATLASGAL contours contain 430 M of material 
The bot tom two panels show the column density (left) and temperature (right) maps fr om t he modified blackbody fitting (Section y 
and |3.4| of the Herschel, SCUBA and ATLASGal maps. The point sources from Section [375] are indicated by x. The contours ind 
column densities of 0.85 • 10 22 , 1.10 • 10 22 , and 1.34 • 10 22 cm and encircle 640, 390, and 220 Mq respectively. The white circle indicates 
the pointing and beam size of the CH 3 C 2 H observations, which indicates a gas temperature of 26 zb 7K. The temperatures from the point 
source fitting are indicated in the last plot. This temperature map is overplotted with isothermal contours of 19 K (black) and 18 K (white). 
These isotherms incorporate 250 Mq and 80 Mq for 19 K and 18 K, respectively. 



1. The global background flux offsets are estimated 
from the large scale Hi- GAL maps and subtracted 
from the maps. 



All data (PACS at 160 /mi, SPIRE at 250 and 
350 /im, SCUBA at 450 /im, and ATLASGAL 
870 /im) were convolved to a common limiting 
beam size (SPIRE 350 /im) and converted to uni- 
form units (Jy/arcsec 2 ). The inclusion of the also 
available SPIRE 500 /im data would force us to 
commonly convolve all other maps to it coarse res- 
olution of around 36". We therefore neglect these 
500 /im data for the temperature mapping. 



3. In the final fitting process an SED was fitted to each 
pixel using a modified black body function, requir- 
ing at minimum three data points for the SED. 

For the final pixel-to-pixel SED fitting a careful prepa- 
ration of the initial data is needed. A complete descrip- 
tion of the data treatm ent and fitting process is g iven in 
lLaunhardt et al. |(|2012| ). The observed globules in |Laun-| 
|hardt et al. | ( |2012| ) are more isolated and suitable back- 
ground flux offsets can be calculated within the EPOS 
maps. Unfortunately, the larger extent of the IRDC, the 
imprint of t he adjacent GMC and infra-red bubble (see 
Section 3.1), and the limited map sizes prevent a rea- 
sonable background determination. To overcome these 
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limitations we utilized Herschel maps from the Hi-GAL 
program which cover a much larger area. We extracted 
the background fluxes from a low-emission region off the 
galactic plane located at 19 h 23 m 26 s .018 +12°37 / 28 / /06 
( J2000) with a box size of 500" . The actual method 
to extract the background fluxes in this fiel d with rela- 
tively low spatial variations is adopted from |Stutz et al.| 
(2009). It is an iterative Gaussian fitting algorithm ap- 
plied to the pixel flux distribution in the selected region. 
The best fit to the background levels is found by omit- 
ting the non- Gaussian components. Finally background 
levels of 0.040 Jy/pixel and 1.434, 1.191Jy/beam were 
c omputed for 160, 250 a nd 350 /im, respectively. 
Battersby et al.| ( [2011 ) used a large scale Gaussian pro- 



file to remove the Galactic background from their ob- 
served Hi-GAL maps toward two different fields before 
the SED fitting. Herschel temperature maps of a sub- 
region c lose to the star - formi ng complex W43 were cre- 
ated by Beuther et al. (2012) without any global back- 
ground subtraction. For lower temperatures (< 25 K) 
both methods just differ by 1 — 2K. Our maps appear 
to be more sensitive to the background determination^] 
however changing our estimated background level by 
±50% does not have a significant impact on the derived 
temperatures. 

For the SED fitting all Herschel maps (PACS 160 /im, 
SPIRE 250 and 350 /im) were convolved to the 



SP IRE 3 50 /im beam, using the kernels from |Aniano et 
aT1 fl2011|) . For the sub-millimeter maps (SCUBA 450 /im 
and APEX 870 jam) a Gaussian kernel was used. For this 
step all maps were calibrated in units of Jy/arcsec 2 . 

In the resulting maps the remaining emission per pixel 
is assumed to be well represented by modified blackbody 



emission: 



S v (v) = n(l-e- T M)B v (v,T A ) 



with 



t(v) = Nh ran R 



(2) 
(3) 

where S v {v) is the observed flux density for a given 
frequency v and solid angle Q (all fluxes are normal- 
ized to Jy/arcsec 2 ), t{v) is the optical depth, B v (y,T&) 
the Planck function, the dust temperature, mn 
the proton mass, R the dust-to-gas mass ratio, and 
N H = 2 • N(R 2 ) + ^( H ) tne total hydrogen column den- 
sity. For the column density map (Figure [5| bottom left) 
we assume that all hydrogen is present as molecular gas 
in G48. For the dust opacity k we adopt the model from 
Ossenkopf fc Henningj (|1994| ) for thin ice mantels and a 
gas density of 10 b cm . For the hydrogen- t o- dust mass 
ratio we adopted a value of R = 110 (e.g., |Sodroski et 



al.|1997 ). To account for helium and heavy elements one 
would cnoose the total gas-to-dust mass ratio 1.36 times 
higher. A modified black body function, as in Equation 
[2] was fitted individually for each pixel using N-r and 
as free parameters with a ^-minimization. The Hi-GAL 
data includes only the 70 and 160 /im maps for the PACS 
bands. To avoid contributions from stochastically heated 

5 We also determined the background offsets from a field with 
low spatial variations and almost no compact emissions inside the 
galactic plane. But the resulting temperature maps from the SED 
fitting were suspiciously high (> 20 K inside the IRDC). Therefore, 
we decided to estimate the background offsets further outside of the 
galactic plane. 



small grains or optic al thick emission, we exclude the 
70 am ba nd (e.g., Battersby et al.|20"TT| Pavlyuchenkov et 
al.|2012[ ). Usually the fluxes in the three remaining Her- 
schel bands (160, 250, and 350 /im) are well constrained 
and provide a good fit to the modified black body func- 
tion. Therefore, every shown pixel in the temperature 
map (Figure [5]) results from an SED fit of at least these 
three Herschelbands. If one of the submillimeter (450 or 
870 /im) fluxes deviate by more than 3 a from the SED- 
fluxes constrained in the other bands it is not used in the 
temperature and column density calculation. 

The temperature map (Figure [51 bottom right) shows 
a steep drop at the boundaries of G48 within the reso- 
lution limits set by the size of the convolution kernels (a 
19 K contour is overplotted on the IRAC map at 8 /im 
and PACS at 100 /im in Figure [3]). In the central region 
of G48, with the strongest sub-millimeter emission, the 
temperatures drop below 18 K (white contour in Figure 
[5] bottom right). The local minima agree with the posi- 
tion of the sub-millimeter peaks PI (17.1K),P2 (16.8 K), 
and EP (17.1 K). In general our observed line-of-sight 
te mperatures are high er than the ones (10-15 K) found 
by |Ormel etH] (|20Q5) for PI, P 2, and EP. However the 
later temperature profiles from |Ormel et al.| ( |2005[ ) are 
obtained fr om fitting observatio nal results t o cor e models 
(similar to |Wilcock et al. 12012 , see Section 4.5). 

To assess the uncertainties 01 our line-of-sight- averaged 
temperatures, we investigated how the change of differ- 
ent parameters influence the resulting temperatures and 
column densities. When altering the subtracted global 
background level, that is subtracted from the Herschel 
maps, by ±50% before the fitting process and using dif- 
ferent beam sizes for the convolution, the resulting tem- 
peratures are found to differ by no more than 0.3 K inside 
the IRDC, where dust emission is detected in the AT- 
LASGAL maps. In the temperature fitting routine, we 
assume mean absolute background flux levels that were 
derived for regions well outside the galactic plane (see 



|Launhardt et al. ~ 2012[ ) and might not be appropriate 
tor this IRDC. For lack of better estimates of the local 
absolute background levels, we have tested the effect of 
simply amplifying those background levels by a factor of 
100. The resulting temperatures and column densities in 
the IRDC were found to change by no more than 0.1K 
and 2.5%, respectively. The effect would be stronger for 
lower column densities, which we do not consider here. 
Hence, we can robustly conclude that the effect of the un- 
certain knowledge in the absolute local background flux 
levels on the derived temperatures and column densities 
in the IRDC is smaller than the uncertainties introduced 
by the calibration, data reduction, and dust opacity mod- 
els, and is thus negligible. 

The total relative uncertainties introduced by calibra- 
tion and data reduction are found to be 2.5% for the tem- 
peratures and 10% for the column densities. The largest 
uncertainties are introduced by the uncertain knowledge 
of the underlying dust opacity model. Therefore, we as- 
sume an overall temperature uncertainty of ±1 K. A more 
detailed discussion of the different contributions to the 



uncertainties can be found in Launhardt et al. (2012). 



It should be mentioned that the temperatures outside 
G48 are not very reliable because of the brightness gradi- 
ent introduced by W51 and N100, as mentioned earlier. 
In the column density and temperature maps (Figure 
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[5| the point source positions and temperatures are indi- 
cated. 

The temperatures derived at the locations of the point 
sources in the large-scale temperature map are dis- 
crepant from the temperatures derived by fitting the in- 
dividual SEDs of said point-sources, the later of which 
were extracted at native resolution. This discrepancy 
is not surprising as it is caused by beam convolution in 
the temperature map and the inclusion of longer wave- 
lengths. 

3.4. Column densities 

To estimate the column densities for G48 we utilized 
three different methods; mid-infrared extinction map 
from IRAC at 8 /im, sub-millimeter dust emission and 
modified blackbody fitting including Herschel and sub- 
millimeter data (see Section 3.3). 
The IRAC 8 /im extinction maps provides us wi th an 



estimate of the optical depth per pixel (see Section 2.2). 
These optical depths were conve rted into column d ensi- 
ties (Figure [5j bottom left) using QVasyunina et al.|2 QQ9): 



JV Ha = 1.086 r/C ext 



(4) 



Where r is the optical depth per pixel and C ext is the 
extinction cross-section. The nature of the mid-infrared 
extinction law below 10 /im in IRDCs and massive star- 
forming regions is still a matter of debate (e.g., |But-| 



iormmg regions is still a matter 01 debate (e.g., |l3ut-| 
ler fc Tan[|2009|). In agreement with findings of |lnde-| 
betouw et al.| ( |2005 ) we used an extinction factor oi 



tiy = 5.5 to parametrize the extinction cr oss-section 
Ce.x.t- We choose C L xt = 4.62 • 10 -23 cm 2 from |Weingart-| 
ner & Draine ( 2001| (case B)). The derived column densi- 
ties also depend on the assumed Ry value and, therefore, 
the used extinction cross-section. The extinction factor 
Rv = 3.1 is more representative of the diffuse ISM and 
would result in significantly higher column densities. An- 
other effect for the 8 /im extinction is the saturation of 
the column density for high values, therefore the masses 
based on IRAC 8 /im discussed further belo w can only be 
lo wer lim its (for a detailed discussion, see |Vasyunina et 
aTJ2009l ). 

The second method to determine the hydrogen column 
densities utilizes the dust emission detected in the 870 /im 
map from ATLASGAL: 



Nn 2 



RF V 



(5) 



where F v is the sub- millimeter flux, R = 110 the gas- 
to-dust ratio, B v (y^Ta) the Planck function for a given 
wavelength (870 /im) and dust temperature (X^ = 17.3 K 
or 19.2 K, see below), mn 2 the mass per hydrogen 
molecule, and ft the beam size. We adopt a mass ab- 
sorption co e fficien t of n = 1.47cm 2 g _1 from [Ossenkopf 



Henning (19941 thin ice mantel, 10 5 cm 3 , as in Sec- 



tion |3^ 

We have three independently derived column density 
maps: IRAC 8 /im, ATLASGAL 870 /im, and the SED 
fitted map combining Herschel and ATLASGAL. Each 
one requires us to make different and independent as- 
sumptions in order to derive column densities; further- 
more, they also vary widely in angular resolution. 

For example, the 8 /im extinction map, which has a 
resolution of ~ 3", requires an assumption of a certain 



extinction factor for the dust extinction cross-section, op- 
tically thin emission and the local background contribu- 
tion (Jbg i n Equation [TJ, but is temperature indepen- 
dent. Furthermore, the residuals of the PSF removal 
for bright IRAC point-sources with compact background 
emission can result in an underestimate of local column 
density (see e.g., source 1 in Figure [5j top left). 

For the ATLASGAL 870 /im map, with a resolution of 
~ 19", always a fixed temperature is assumed, it requires 
optically thin emission and a dust opacity choice at the 
wavelength of interest. 

While the SED-fitted map also requires a dust model 
assumption, it contains information from multiple wave- 
lengths and simultaneously derives 7Vh 2 and temperature 
for each pixel, under the assumption that the mass is 
well-represented by a modified black body and assum- 
ing a fixed spectral index j3. Therefore the derived col- 
umn densities can be considered more robust compared 
to the two previous methods. However, the SED-fitted 
map also has the lowest resolution of the three methods 
considered here, with final convolved beam size of 25", 
equivalent to ~ 0.3 pc at the distance of G48. Here we 
discuss the integrated masses for these three independent 
methods. 

The total mass of G48 can be calculated within a con- 
tour of interest using the column densities mask (shown 
in Figure [5| : 



m H2 Ni(K 2 ) 



(6) 



where a px is the pixel area, mn 2 and A^(H2) the column 
density per pixel. 

The masking contours for the IRAC 8 /im extinction 
map were chosen by visual inspection. This contour re- 
sults in a threshold of 32.15 MJy/sr, which correspond to 
a column density of 1.34 • 10 22 cm -2 . A total H 2 mass of 
450 M is enclosed in these contours (see Figure [5| top 
left panel). 

We assume a temperature of 17. 3 K (the median point 
source temperature, see Section |3.5|) for the ATLAS- 
GAL derived column density map (Figure [5] top right 
panel) . The column density threshold is adopted to be at 
3xthe RMS noise level (~ 50mJy/beam), corresponding 
to AT(H 2 ) = 0.85 • 10 22 cm" 2 . The resulting total mass 
inside this column density contour is M(H 2 ) = 430 M . 
For comparison, the median Herschel derived tempera- 
ture (Figure [5] bottom right panel) in this same region is 
19.2 K. If we adopt this temperature instead of 17.3 K, 
as assumed above, the column density threshold becomes 
N(R 2 ) = 0.75 • 10 22 cm" 2 and the total enclosed ATLAS- 
GAL mass is M(H 2 ) = 675 M . 

For reference, applying a column density threshold 
of 1.34- 10 22 cm -2 (derived from the from IRAC 8 /im 
column density map) is equivalent to applying a 18 K 
isotherm to the Herschel derived column density map, 
and results in thresholds that encompass approximately 
the same spat i al ext ent of the three SCUBA peaks from 
lOrmel et all (120051). The mass enclosed within these 
thresholds differ by only 20 — 30% compared to the total 
mass (320 M ) of these three peaks. 

The total mass for G48 derived from the 8 jam extinc- 
tion (450 M ) and the 870 /im (430 M ) emission are re- 
markably similar. Furthermore, the average column den- 
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sity threshold for both maps corresponds to a Herschel 
derived total mass of M(H 2 ) = 390 M , also in excellent 
agreement with the other estimates. Considering the un- 
certainties in the various assumptions, such as the dust 
properties, uncertainties in background flux levels, and 
variations in beam sizes we find that our mass agreement 
for G48 is strikingly good across different methods and 
compared to previous work. 

In general the column density for the diffuse emissions 
in G48 are comparable for all three methods. For EP, PI, 
P2, P3 and P4 the peak column densities agree within 
18% for the ATLASGAL and the Herschel maps. How- 
ever, the peak column densities derived from the 8 /im 
extinction are significantly higher then the fitted Her- 
schel maps. The peak columns can be higher up to 75%. 
This might be an effect introduced by the higher spatial 
resolution of IRAC. 

3.5. Point source Classification 

In Table [2] we present the fluxes of the 13 point sources 
which were detected in all three PACS bands. Except for 
source 8, all PACS detections have count erparts in MIPS 
24 /im. Source 1 3 was not included in the |van der W iel & 
Shipman| ( |2008| ) sample, but is detected as a point source 

If one allowed point source de- 
bands this would add only one 
This source is also detected 



at 24/xm (E21h Figure m 
tection in only two PAC 
more obvious point source, 
at 24 micron (El) but appears as a diffuse blob of emis- 
sion at 100 micron, not well-matched to the characteris- 
tic PSF we expect. However, the projected distance to 
Source 13 (E2) is ~ 0.24 pc which is close to the median 
separation discussed in Section [4j and therefore would 
not change the outcome of this discussion. All detected 
sources are aligned along the mid-infrared filaments of 
the IRDC. Even though the emission gradient towards 
the GMC (on the eastern edge of our maps) is brighter 
and contains some compact emission features, no PACS 
point sources a re de tected outside the IRDC filaments. 

As in Section [33} we used a modified blackbody func- 
tion to constrain the parameters of th ese point sources 
(Figure |6j ), utilizing the dust model by Oss enkopf fc Hen^| 
|ning| (|1994[ thin ice mantles, 10 5 cm -3 ). The sources de- 
tected at 24 jam cannot be fitted by a single component 
fit, since the warmer part of the envelope in the direct 
vicinity of the evolving protostar is contributing signifi- 
cantly to the 24 /im fluxes. Due to the rapidly increas- 
ing optical depth towards the innermost regions around 
the point-sources we only observe surface contributions 
of surrounding layers of warm dust. The calculations of 
mass and temperature assume optically thin material. 

The 24 /im data are hig hly sensitive to geometric effects 
in protostellar cores (see |Dunham et al.| |2008; Rag an et] 
al. 2012[ ). A proto-star m which outflows have carved 



a cavity in its natal core will be detectable at 24 /im if 
our line of sight is aligned with the cavity. In partic- 
ular the shocked molecular hydrogen, visible as green 
structures in Figure [l] (right, bottom panel), traces high 
velocity outflows. They originate from the rotational S9- 
Sll transition of excited H 2 in the IRAC 4.5 /im band. 
They are generally referred to as either "green and fuzzy 
features" (Chambers et al. 200 9]) or "exte nded green ob- 
jects" (EGOs, [Cyganowski et al. |2008[ ). Because the 
24 /im flux arises primarily from the warm material im- 
mediately surrounding the proto-star (and not the cold 



outer region of the core), we do not include the 24 /im 
flux in our SED fit to the PACS data but instead fit 
a second modified blackbody component (see Figure [6]). 
The best fit to the SED results from the sum to these 
two components. For Source 8, we used one modified 
blackbody function. The fit parameters are presented in 
Table 12 

The four most massive sources in our sample are lo- 
cated in the central part of the IRDC in the direct vicin- 
ity of the sub- millimeter peaks PI and P2. They cover 
a mass range from 5.4 to 14.3 M . From these sources, 
Source 8, is the only PACS source without a MIPS 24 jam 
counterpart. Source 1 and 2 are the brightest and most 
massive sources in the sample. They are associated with 
the SCUBA peak PI and have a close separation of 
~ 12". Source 1 appears on top of a compact and diffuse 
background emission complex. Only a faint source with 



0.9 Mq is as sociated with the EP-SCUBA-peak of Ormel 



et al. (2005). The point sources associated to P3 in the 
southern part of G48 have only ~ 1 — 2 M . 

In Figure [7] we present the bolometric luminosities 
(IRAC to PACS), the envelope masses (< 0.15 pc), and 
the temperatures derived from cold components of the 
modified blackbody fits. In the luminosity-mass plot the 
dotted line follows the relation M env ~ and corre- 
sponds approximately to a threshold, where 50% of the 
initi al core mass being converted into stellar mass 
(c.f. |Andre et al.||2000| ; |Bontemps et al.||1996[ ). Sources 
with M env > M* and a high sub-millimeter to bolomet- 
ric luminosities rati o I/>35o //m / Lpoi > 0.005 are consid- 
ered to be Class ( | Andre et al. 2000). Source 1 and 2 
as well as Source 4 and 8 are detected as single SPIRE 
sources. But the fitted blackbodies indicate that the sub- 
millimeter fluxes should be equal for both sources in each 
case. For each source we calculate the sub- millimeter 
luminosity £a>350/i™ with 50% of the corresponding 
SPIRE fluxes. All four sources mentioned before are 
good Class candidates since L>s50fim/Lb o i > 0.005 and 
M env > M*. Source 1 and 2 have the highest enve- 
lope masses, an d, as indicated by evolutionary tracks 



(adopted from | Andre et al. 2008| ), both could evolve 
into high-mass Zero Age Main Sequence stars (ZAMS) 
with Mzams ^ 8M . These evolutionary tracks as- 
sume exp onentially dec l ining accretion rates. A recent 
study by Davies et al. (2011), comparing observations 
to various accretion models, showed that models with 
increasing accretio n are in best agreemen t with the ob- 
servational results. |Molinari et al. (2008) used a model 
with increasing accretion rates to create similar evolu- 
tionary tracks. According to these tracks Source 1 and 2 
would only evolve into intermediate-mass stars (< 8M ). 
However, it is still unclear on which spatial range a proto- 
star is accreting material. Source 1 and 2 are associated 
to the sub-millimeter peak PI. Due to the resolution 
of SCUBA compared to PACS a larger area is traced. 
Given the li ght-to-mass ratio for PI derived by Ormel 
et al. (2005), still two massi ve stars (Mzams ^, 8 M fi ) 



could form m the model from Molinari et al. (2008). As- 
suming these sources are in a very young evolutionary 
stage, they are completely dominated by the accretion 



2009). 



lumin osity (Krumholz et al. 2007 Hosokawa & Omukai 
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Fig. 6. — The modified blackbody fits for the detected sources (see Section [375]) . 



TABLE 2 
Point sources 



Source 


R.A. (2000) 
[hh:mm:ss[ 


Dec (2000) 
[dd:mm:ss] 


mJy 


F70 ixm 

mJy 


mJy 


mJy 


T L 
K 


L L 

L 


M L 

M 


L H 

L 


Comments^) 


1 


19:21:49.65 


+13:49:31.26 


102.7 


1764 


4086 


6931 


18.3 


48.2 


14.3 


36.1 


S5, PI 


2 


19:21:50.24 


+13:49:39.28 


241.1 


3702 


6931 


10712 


17.9 


35.5 


11.9 


38.0 


S6, PI 


3 


19:21:50.26 


+13:48:18.96 


119.5 


1730 


3659 


4745 


20.9 


17.7 


2.4 


14.6 


S2, P3 


4 


19:21:47.95 


+13:49:20.37 


2.8 


460 


1807 


4717 


17.2 


15.1 


6.3 


2.8 


S15, P2 


5 


19:21:50.13 


+13:47:57.00 


27.7 


474 


902 


1474 


19.8 


5.5 


1.0 


4.5 


SI 


6 


19:21:40.76 


+13:51:17.42 


102.7 


365 


494 


693 


16.8 


2.4 


1.2 


6.5 




7 


19:21:45.58 


+13:49:20.51 


8.7 


299 


842 


1370 


19.5 


4.8 


1.0 


1.7 


S12, EP 


8 


19:21:48.48 


+13:49:25.26 




210 


865 


4668 


17.4 


13.3 


5.3 




P2 


9 


19:21:53.01 


+13:48:19.15 


34.7 


139 


174 


357 


16.5 


1.1 


0.6 


2.8 


S3 


10 


19:21:36.30 


+13:52:20.23 


31.8 


132 


780 


2209 


16.8 


7.9 


3.8 


1.6 


Sll 


11 


19:21:35.95 


+13:51:59.97 


115.5 


121 


177 


521 


15.6 


1.7 


1.3 


4.4 


S9 


12 


19:21:40.73 


+13:50:54.91 


20.2 


115 


303 


894 


15.6 


3.2 


2.4 


1.5 


S13 


13 


19:21:31.93 


+13:52:51.82 


4.6 


61 


131 


516 


14.5 


1.7 


1.9 


1.0 




Note. — (a) Sx are 
millimeter peaks. 


the associated MIPS 24 /xm 


from (van 


der Wiel 


& Shipman| <|200£ 


|, Px 


and EP 


are the associated sub- 



3.6. Molecular Line Emission 

To determine the parameters from the mapping line 
data, we used the GILDAS/CLASS package. As men- 
tioned before we have line detections (CL> 99.73%) as- 
sociated to the IRDC in the HCO+, N 2 H + , HCN, HNC, 
13 CO and C 18 maps . While the line fitting for HCO+, 
HNC and CO is done straightforward with a single Gaus- 
sian profile, N2H + and HCN have a resolved hyperfine 
structure with multiple components. For N2H + a model 
profile containing seven hyperfine- structure components 
was fitted utilizing the GILDAS/HSF routine. The pre- 
sented N2H + maps containing the integrated line flux, 
line width and velocity correspond to the line center 



6 http://www.iram.fr/IRAMFR/GILDAS/ 



of the strongest component (93.17378 GHz , J = 1 0, 
Fi = 2 — >► 1, F = 3 ->> 2). The hyperfine-structure of 
HCN shows a triplet of components, and in contrast to 
N 2 H + , the relative li ne strength of each component varies 
for di fferent IRDCs ( |Afonso et al.||1998| Vasyunina et al. 
2011). Non-L TE effects especially in tne 1-0 lines seem 
common (see Loughnane et al. 2012). Given the noise 
level of the HC1N map we could not retrieve an averaged 
profile. The fitting of the physical parameters is done on 
the strongest component. 

The derived integrated line fluxes, velocities and line 
widths are compiled in Table [5] This Table gives also 
the effective radius (see belowj, the total column den- 
sity, the abundance and the traced mass. HCO+ and 
HCN are good tracers for dense gas because of their high 
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1x10" 1 1x10° 1x10 1 1x10 2 1x10 3 1x10 4 
L b0 , Hi 

Fig. 7. — The core temperatures (upper panel) and cold envelope 
mass (lower panel), derived from the SED fitting of the PACS point 
sources (< 0.15 pc), are plotted against bolometric luminosities. 
Bolometric luminosities are derived from IRAC, MIPS and PACS 
bands. Th e numbers identify the point sources characterized in 
section [33] The light-to-mass ratio derived from the sub-milliliter 
observation for different peaks (P1,P2 and EP) is plotted as well. 
In this observations the area of the traced mass reservoir is larger. 
The dashed lines represent the evolutionary tracks for YSO from 
|Andre et ah] (2008). The finale stellar masses (ZAMS) of those 
tracks are given above the lower axis. 



TABLE 3 

Clumpfinder parameters 



Molecule 


lower threshold ( a ) 


step ( a ) 




% 


% 


N 2 H+ 


5.0 


5.0 


HCO+ 


5.8 


5.8 


HCN 


8.5 


3.4 


HNC 


18.6 


9.3 


13 CO (1-0) 


35.2 


2.3 


13 CO (2-1) 


31.4 


6.3 


C 18 Q (2-1) 


8.8 


5.5 



Note. — ( a ) The minimum detec- 
tion thresholds and the contrast steps for 
clumpfinder are give in relation to the peak 
value in the corresponding map. 



critical densities. N 2 H + has a lower critical density (see 
Table [I]) , but is not strongly effected by depletion from 
freeze-outs on grain surfaces. For CO the critical den- 
sity is, depending on the transition, one or two orders 
of magnitudes lower than for N 2 H + . The total column 
density for a certain molecular species can be calculated 



TABLE 4 

Upper limits for non detected line species for 
MOPRA 



Line 


line integral intensity ( a ) 


column density 




[kms _1 K] 


[lO^cm- 1 ] 


CH 3 C 2 H 


0.50 




H 13 CN 


0.53 


3.8 


H 13 CO+ 


0.52 


0.7 


SiO 


0.51 


1.2 


C 2 H 


0.52 




HNCO 


0.48 


6.4 


HCCCN 


0.48 


2.0 


CH 3 CN 


0.50 


2.9 


13 cs 


0.53 


5.9 



a Upper limits are the integrated intensity a Gaussian- 
shaped line with a FWHM of lkms -1 would have at a 
peak intensity level of three times the RMS noise. 



with 



tot 



%TtA ul kv 2 

N u Qrot{T) 



T m bdv 



exp(-E u /T) 



(7) 



(8) 



where A u \ is the Einstein coefficient, v the rest frequency 
of the transition, N u and E u the upper state column 
density and energy, Q ro t the rotational partition func- 
tion and g u the statistical weight of the upper level. The 
general physical parameter^Jfor the observed lines are 
given in Table [TJ The abundance for a molecule is given 
in relation to the H2 colum n de nsity derived from the 
SED-fitted map from Section 3.3 The temperature esti- 
mates needed to calculate the different molecular abun- 
dances were obtained from the same SED-fitting. The 
line emissions from HCO + , HCN, HNC and 13 CO can be 
optical thick, therefore the column densities and abun- 
dances provide only lower limits. 

For the molecular line maps without any 3a detection 
the upper limits for line intensity and column density, 
based on th e noise estimates, ar e given in Table |4j 

Following |Ragan et al.| ( |2Q06|> the mass of each clump 
observed in the line maps is estimated using an assumed 
distance of 2.6 kpc, an approximate size based on the ex- 
tent of molecular emission, the molecular column densi- 
ties, and approximate abundance calculated at the peak 
of absorption. 

In Figure [4] the SPIRE 250 /im images are shown and 
overplotted by the different line intensity maps (HCO + , 
HCN, HNC an d CO). We used the clumpf ind 



NH+, 



(Williams et al. 1994) package to identify the different 
clumps in the molecular tracers. Four different clumps 
were found by clumpfinder, hereafter named C1-C4. The 
detection parameters for clumpfinder are given in Table 
[U CI, C2 and C3 are detected in HCO+, HCN and 
N 2 H + . An additional clump C4 is detected in HCN. In 
the 13 CO (1-0) GRS map only CI is detected. While CI 
and C2 are detected in the HERA CO (2-1) maps and in 
HNC. Because clumpf ind does not make any assumption 

7 The parameters can be obtained from the Leiden Atomic and 
Molecular Database (LAMBDA). 
|http://www.str w. leidenuniv.nl/~moldata/ 
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about the clump morphology an effective circular radius 



R, 



eff 



(9) 



is given in Table [5j where A is the area fitted by 
clumpf ind. Because the resolution of the HERA obser- 
vations is three times better than MOPRA, the emission 
is less diluted and we therefore derive smaller R e ff- All 
line maps show the strongest emission in the central part 
of the IRDC. The maps HCO+, N 2 H + , HNC, and 13 CO 
show an elongated structure in east- western orientation. 
HCN shows two distinct peaks (CI, C4); these peaks 
are shifted by ~ 30" south-west of the compact SPIRE 
emission. All other MOPRA molecular tracers peak at 
CI, and while the maps appear shifted from the com- 
pact SPIRE emission, it is within the pointing error of 
MOPRA. The peak positions for the HERA CO maps 
are much better aligned with the SPIRE peaks. In con- 
trast to our observed HCO + (1-0) emissi on, the HCO + 



(3-2) observation from Ormel et al. (2005) show multiple 
emission peaks. Their eastern emission peak is centered 
on the SCUBA clump, but the western emission is more 
diffuse. With the given spatial resolution of MOPRA 
(38") it is likely that we only observe a single integrated 
emission peak for HCO + from these multiple unresolved 
emitting sources. 

About 1.'3 southwards of the central region another 
fainter emission peak (C2) can be found for HCO + , HCN, 
HNC, N 2 H + , 13 CO(2-l) and C 18 0. It is not coincident 
with the compact SPIRE and sub-millimeter emission in 
this region. In the integrated line flux map for 13 CO(l- 
0), no compact emission clump was found by clumpf ind. 
The emission appears just as a diffuse extended lobe from 
the central peak. In the north-western part of the IRDC, 
where only diffuse emission appears in the SPIRE and 
sub-millimeter regime, the line data are much more noisy. 
For HCO + , HCN and N 2 H + a compact emission peak 
(C3) can be identified, while 13 CO shows a elongated 
diffuse emission lobe. 



The three SCUBA clumps modeled by |Ormel et al. 
( 2005| ) incorporate a mass reservoir of ~ 320 M . The 
masses derived for this region from the HCO + , HCN, 
HNC, N 2 H+ and 13 CO observations differ by up to 20%. 
Given the uncertainties for the mass estimates, this is a 
good agreement. 

CH 3 C 2 H is a go od temperature prob e in dense molec- 

We detect four 
= 0,1,2,3 with 



ular clouds (e.g. 



Bergin et al. 
7 



1994) 

8), 



K 



CH 3 C 2 H transitions! J = (9 
IRAM single pointing observations positioned on the sub 
millimeter peak PI (as indicated in the bottom left panel 
of Figure pi . From the measured integral line intensi- 
ties we derived the column densitie s for every pa r ticula r 
transition using formulae 1-2 from iBergin "etaD (fT994l). 
Finally we computed the excitation diagram and esti- 
mated the excitation temperature and column density 
in the region. The resulting temperature is 26 ± 7 K and 
the column density (3.0 ± 1.5) • 10 12 cm -2 . All necessary 
values, like E u - the upper state energy, g u - degen eracy 
and others, were taken from the CDMS data base ( |Muller| 
[eTaL|[200T| ) 

In Figure [8] (bottom panel) the position- velocity (p-v) 
plots are shown for 13 CO and C 18 (2-1). In the p-v plot 
for C 18 only a single velocity component is present. For 




Fig. 8. — In the top panel the contours of the line intensities 
for 13 CO (left) and C rs O (right) is plotted over the inversed 
PACS 70 /xm image. For 13 CO (HERA) the blue contour is in- 
tegrated between 30.0 — 34.0 kms -1 and the red contour between 
36.0 - 40.0 kms -1 . The integration in the C ls O HERA map is 
done over a more narrow range, due to higher noise level, with 
32.0 - 33.5 kms -1 (blue) and 34.5 - 36.0 kms -1 (red). 
The middle plot shows the linewidth for the HERA CO data. The 
overplotted contours are the isotherms of 17.5, 18.5, 19.0, and 
19.5 K (see Figure [5}. 

The bottom plot shows the position- velocity diagram for 13 CO and 
C ls O. The spatial orientation of the p-v-plot is indicated in the top 
panels by a black vertical line. 



13 CO two components are present at a 50" offset, which 
refers to the spatial position of the Source 1 in the PACS 
points- source catalog and the sub- millimeter peak PI. 
Due to the small spread of just ~ 1.5 kms -1 it is ques- 
tionable that this structure originates from an outflow 
cone. Other possible explanations could be a large rota- 
tional envelope or just two core-like dense CO structures 
along the line-of-sight. With the lack of additional ob- 
servations with higher sensitivity and spatial resolution, 
no robust conclusion can be drawn. The linewidth maps 
in Figure [8] are fitted with a single Gaussian. A single 
Gaussian profile is also fitted there two components are 
found in the 13 CO map (Figure [5J middle plot). The 
resulting does not differ more then 0.5 kms -1 from a fit 
only on the strongest component. 

With typical isotopic abundance ratios of 
12 C1/ [ 13 C1 =60 and [ 16 01/[ 18 Q1 = 500 (|Wilson 



Rood 



1994 



Lodders et al. 2009) the HERA l6 C0 and 



C O abundances correspond to 12 C 16 abundances of 
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~ 1 — 1.5 • 10 -5 . The expected CO abundance at a given 
Galactocentric distance ( Dqc) can be compu ted from 
the following relationship ( |Fontani et al. 2012 and ref. 
within): 

Xg 16Q = 8.5 • HT 5 exp(1.105 - 0.13L> G c/£>k pc ) (10) 

For G48 with Dqq = 7.1 kpc this results in an expected 
CO abundance of X^ 16Q = 1.0 • 10 -5 thus consequently 
Xg 18Q = 2.0 • 10- 7 and X? 3co = 1-7 ■ 10~ 6 . Comparing 
this expected values with the observed CO abundances 
gives the depletion factor: 



f D = 



Y E 
A CO 

*8 b o s 



(11) 



Using the preceding isotopic abundance ratios results in 

3 CO o o ___ _1 ^C 18 _ ^ rj 



c 18 o 
3.3 and f D 



average depletion factors of f D 
for G48. 

We detected C0 2 ice in the Spitzer spectra. In Fig- 
ure [2] the Spitzer /IRS spectrum for Source 2 is shown. 
The observed shape of the C0 2 ice absorption feature 
at 15.2 fim is different for the source spectrum compared 
to the background spectra; the on-sources profile is sig- 
nificantly broader and different in shape compared to 
the nodded background profiles. This implies differ- 
ent ice temperatur es and mixing ratios (with e.g., H 2 
and CH3OH, see [White et al.||2009| ). The low spec- 
tral resolution of the 1RS/LL2 spectral order does not 
allow a proper fitting of temperature and mixing ratio 
to the laboratory data. Most likely the source C0 2 pro- 
file originates from ice coated grains in the dust enve- 
lope of the embedded proto-star. The presence of C0 2 - 
ice does not provide direct evidence for CO freeze-out, 
since the conditions under which the surface conversion 
9 + H takes p lace are still under de- 
Garrod fc Pauly|2011 ). In the gas-phase C0 2 



CO + OH CO 



bate (e.g., 

is predominantly produced by the following reaction 



HCO + O 
CO + OH 
CO + O 



C0 2 
C0 2 
C0 9 



H 
H 
H 



To the last two reactions an activation energy barrier of 
80 K or higher is assigned (Garrod & Pauly 2011 ). There- 
fore, gaseous C0 2 creation becomes effective in cores 
which already have an internal heating source. The C0 2 
ice cannot be easily hydrogenated on dust grain surfaces 
though, so it will remain a terminal species of the CO 
evolution in lukewarm regions (Dmitry Semenov 2012, 
private communication) . 

4. DISCUSSION 

G48 is a quite isolated IRDC at a distance of 2.6 kpc. 
Consequently, G48 provides nonpareil observations of the 
collapse and fragmentation of a IRDC without the effects 
of external feedback, except for the diffuse ISRF. 

We have derived line of sight averaged Herschel col- 
umn density and temperature maps for this source, and 
find minimum temperature of 16.8 K, and outer temper- 
atures of ~ 23 K. The corresponding column densities 
range from ~ 3 • 10 22 cm -2 at in central part of G48 to 
~ 1 • 10 22 cm -2 at the outer edges where our sensitivity 
drop sharply. 



In this section we will discuss these results in combi- 
nation with a wealth of other observations with the goal 
placing this unique source in an evolutionary context and 
determining its physical properties. 

4.1. Molecular line data, CO freeze out 

The comparison of the line data for G48.66 with 
other IRDCs show results for N 2 H + which are similar 
to other IRDCs, while the HCO + results are different. 
The average line width of HCO + is 2.5 kms -1 which is 
lower compared to other IRDCs, with typical values of 



3.2 — 3.3 kms (Vasyunina et al. 



2011 



Zhang et al. 



2011). For N 2 H + we observed an average line width of 
1.8 kms -1 . Similar mean line widths were observed for 



other IRDCs with 1.8 kms" 1 byJVasyunina et al.| ( |2011 ) 
and 2.0 kms -1 by Ragan et al. (2006). The mean col- 
umn density of 2.4 • 10 12 cm -2 for N 2 H + in our sample 
is similar to the column density found for other IRDCs 
by Gibson et al. (2009). Comparing the amount of CO 
freeze out one finds that the depletion factor is strongly 
dependent on the chemical age of the IRDC, h ence the 
depletion varies drast ically for different IRDC s. |Hernan-] 
dez et al. (2011) and Hernandez et al. (2012) reports a 
depletion factor of 5 for G035.39.-0033 which is a young 



Miet- 



IRDC of > 7 : 10 4 yr similar to G48. In contrast 
tinen et al. (2011) found no depletion for their sample 
and argued that this it is connected to a more evolved 
stage of their IRDCs where the increased star forma- 
tion activ ity leads to sublimat ion of CO from the grain 
surfaces. Fontani et al. (2012) reports higher values of 
fr> = 5 — 74 lor several southern IRDCs, the reason for 
this higher depletion factors is not completely clear, but 
might be explained by tracing denser regions of the cores 
by utilizing higher C 18 (3-2) transitions, compared to 
our C 18 (2-1) observations. In the gaseous inter clump 
is destructed i n the following reaction : 
> HCO+ + N 9 ( |Bergin fc Tafalla||2007[ ). 



medium N 2 H + 



N 2 H+ 



CO 



The presence of N 2 H + is consistent with the low CO 
abundance induced by the freeze-out on the dust-grain 
surfaces. 

By c omparing their HC O + observations to clump- 
models Ormel et al. (2005) found an increase of turbu- 
lent motions towards the inner part of the clump. Their 
observations have higher resolution by a factor of two. 
Therefore, our observations are not good enough to trace 
down turbulence in the same way, but the increasing MO- 
PR A line widths in wards of the c l umps are in agreement 
with the findings of |Ormel et ah] (|2005 ) . Our HERA CO 
data have a higher resolution than their observations, but 
due to the low critical density of CO we trace a larger col- 
umn that might be dominated by envelope material also 
along the line of sight. Therefore, CO does not qualify 
as a good turbulence tracer in dense environments as in 
our observations. 

4.2. Cloud stability 

To assess the dynamical state of the molecular clumps 
we calculate the virial mass, using the line width Av: 

M virial = k 2 RAv 2 (12) 

where as fc 2 = 126 for a n ~ 1/r 2 profile (|MacLaren et al. 



|1988[ ), R is the effective radius of the molecular clumps. 
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and Av is the FWHM of the respective line in kins" 1 . 
We also calculated the virial parameter a which is given 
by the ratio between the kinetically energy E^ n and the 
potential energy E pot : 



line mass significantly (Fiege & Pudritz 2000). How- 
ever, the higher line mass ol 385 lor the more evolved 



2£^ki] 
E, 



M v 



pot 



M 



(13) 



A value of a « 1 represents a virial equilibrium, while 
1 < a < 2 is present for pressure confined clumps 
(Bertoldi & McKee 1992), a > 2 indicates a gravita- 
tionally unbound clump, and for a <C 1 the clump is 
collapsing. The virial mass and the virial parameter are 
presented in Table |5j We see a large range, from less 
than 1 to over 7, in virial parameters depending on the 
species, making it difficult to assess the stability of each 
clump with certainty. However, while CI seems to be 
off equilibrium and just confined by the surface pressure 
(aci = 1.9 =b 0.7), C2 and C3 appear to be gravitation- 
ally unbound (ac2 = 3.6 ± 1.3 and ac3 = 3.6 ± 0.4, 
respectively) . 

Another crucial parameter is the virial line mass. An 
isothermal, nonmagnetic, self-gravitating filament in the 
absence of external pressure is in hydrodynamical equi- 
librium as long as the mass per unit length does not 
exceed m v i r = 2a 2 /G, where a is the velocity dispersion 
in the gas. If the filament is only thermally supported 
the veloc ity dispersion is given by the thermal sound 
speed c s ( Ostriker||1964 ), while in the case of additional 
turbulence supp ort the velocity dispe rsion is given by 
a = Av/V8 In 2 ( [Fiege fc Pudritz|2000[) , where Av is the 
velocity dispersion. The critical line mass for a thermally 
supercritical filament in G48 is M^f m = 24M /pc, 
whereas the critical line mass calculated by the mean 
linewidth of G48, to account also for turbulent motions, 
is = 71M /pc. To highlight the derived linemass, 

we straightened the fil ament of G48 utilizi ng a set of non- 
uniform cubic splines ( Kocsis et al. |1991 ) and calculated 
the mass in each pixel. Binning the pixels results in an 
average line mass of 174M /pc (see Figure [9]). The ob- 
served line masses are ~ 2.5 larger than the critical line 
mass, indicating that the G48 filament does not contain 
enough internal support to counter the collapse while 
the contributions from magnetic fields are expected to 
be negligible (see below). Line masses observed for other 
IRDC seem to exceed M^^ m . However, for strongly 
fragmented cores with significantly lower mass in Tau- 
rus line masses similar to M^ e J m : hence predominatel y 
thermal support, were observed (Schmalzl et al. 2010). 
In contrast to G48 several IRDCs show line masses be- 
low or similar to . This cannot be explained by the 
evolutionary stage of the IRDCs alone. The young and 
starless Coalsack cloud and the adjacent G3 04. 74+01. 32 
show line masses of less or equal to (| Beuther et | 



— jj — | — J. j a - - - - ime 

al. 2011 Miettinen 2012 respectively). Also Hernandez 
20121 found " 



eT^PPBl ) found the young H filament ol GU35.39- 
0.33 to be virial undercritical, although these results rely 
several uncertain morphological assumptions about the 
IRDC. G48 and the Nessie IRDC are more evolved, show- 
ing signs of more active star formation such as "green 
and fuzzy" features, but the latter region still shows line 
widths a factor of 5 below the criti cal value for turbu- 
lence support ( | Jackson et al. 2010). Furthermore the 
presence of magnetic fields does not change the critical 



Orion A cloud ( Jackson et al.|20f0 ) supports the view of 
increasing line masses along the evolutionary sequence. 

4.3. Fragmentation and point source separation 

In the central part of the IRDC a concentration of 
point sources is found. The projected median separa- 
tion of PACS-point-sources is ~ 0.3 pc, which is approxi- 
mately the average separation in the central part of G48. 
This separation is about three times smaller than in other 
IRDC^] The median separation in the Rl fragment of 
the Coalsack dark cloud is even an order of magnitude 
higher compared to G48 ( Beuther et al.||2011 ), although 
here this study uses gas/ dust peak separations instead 
of point source separation. The key question is how to 
explain these different values. One could assume that we 
just suffer from a projection effect in G48. Considering 
we see an elongated filament and given the observed mor- 
phology it is a feasible assumption that G48 is aligned 
approximately face-on towards the line-of-sight. Even if 
one would assume that G48 is aligned off by 7r/4 this 
would just account for a factor of y/2. In contrast to 
G48 the mass reservoir in the IRDCs mentioned above is 
significantly higher; these other dark clouds do not share 
the isolated nature of G48. While external feedback and 
the overall evolutionary stage of the parent cloud might 
have an important impact on the fragmentation, evalua- 
tion of the entire h igh-mass part of the EPoS sample in 



Ragan et al. ( 2012[ ) show that the average point source 



separation is not correlated with the cloud mass (Pitann 
et. al, in prep.). Recent molecular line observa tion of 
G28.34+0.06 found a separation of just 0.16 pc ( |Wang| 
|et al.||2011[ ) to 0.19 pc ([Zhang et a l. 2009). The interler- 
ometric data analyzed in Wang et aT" p011| ) show sub- 
structures on a ~ 0.01 pc scale inside the cores. 

Given the wide spread of source separation it should 
be discussed whether G48 exhibits the same mode of col- 
lapse as the IRDCs with larger separation or if some sub- 
stantially different physical processes can account for the 
smaller separation. To address this question we compare 
the observational results to analytical considerations. A 
first simple analysis is to compare the source separation 
to the Jeans length. The separation in G48 is still several 
times larger than the thermal Jeans length (Aj ~ 0.08 pc 
at 10 5 cm -3 , 19 K). Therefore, a thermal Jeans instabil- 
ity is not the dominant fragmentation mode. Fragmenta- 
tion below this spatial scale has been observed for several 
mass ive cores (e.g., Hennemann et al.||2009 Wang et al. 



2011) 



A way to assess the spatial separation scale is to con- 
sider the collapse of a cylinder (with a radius R c ) of an 
incompressible fluid. It is predicted that several cores 
form, period ically placed at a characteristi c length scale 
Afrag = Rc ( |Chandrasekhar fc Fermi| [l953 ) . For an infi- 
nite cylinder ol isothermal sell-gravitating gas the char- 
acteristic fragmentation scale is Af rag = 22H, where H is 
the isothermal scale height H = c s j '(47r /Gpc) 1 ^ 2 and c s 
the thermal sound speed, G the gravitational constant 
and p c the central density in the cylinder (|Nagasawa| 



The separation i s about 0.9 pc for t he IRDCs G011.ll -0.12 
Henning et~aT||2010t an d G049.40-0.01 pCang et al.||20"TT] >; and 
1.75 pc for G3U4.V4-1.32 JMiettinen|2012) . 
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Mass per pixel [MgJ 




ima in their temperature maps are aligned with the col- 
umn density peaks as in G48 and the subsequent studies. 
The minimum temperatures for G030.49+0.96 are lower 
compar ed to G48 and t he tw o fol lowing studies . The 
work of |Battersby et al.| ( |20TT1 ) and |Beuther et al J ( [2012]) 
show similar t empe rature ranges as lor G 48 . W hile | B at- 1 



Fig. 9. — We present the line masses of the straightened filament 
of G48. The upper panel shows the mass per pixel. The lower plot 
presents the masses summed over the width of the filament (hence 
the column of the upper plot). In the middle plot the masses are 
binned to 1 pc. The average line mass for G48 is given by the 
dashed line. 



1987). When a isothermal cylinder with a given radius 



is embedded into a uniform medium the separation is 
Af rag ~ 22 H for R ^> H , but for th e opposite case R <C H 
the length scale is Af rag « 11 R ( [Jackson et al.|[2 QlQ). 
Given the measured temperatures the typical thermal 
sound speed is c s = 0.23 kms" 1 in G48. With 0.5- 0.8 pc 
the width of the filament is much larger than the scale 
height. The calculated length scale Af rag = 22H w 0.3 pc 
(with n = 10 5 cm -3 , c s = 0.23 kms -1 ) is the same as the 
measured separation in G48. However, Af rag is affected 
by uncertainties, since the real value of the central den- 
sity is unknown. Moreover, the most regular positioning 
of sources occur in the central part of the IRDC along a 
clearly changing temperature profile which is a distinct 
derivation from the isothermal assumption. 

4.4. Filament width 
The t ypical filament width seen in Figure [9] is 0.5 



(2011 ) find typical filament widths of 0.08 — 0.1 pc in two 
different low-mass star-forming regions in the Gould-Belt 
Herschel GTO project. These widths are a factor ~ 3 
smaller compared to the SPIRE 350 /im. IRAC provides 
a much better resolution, but using the IRAC column 
density contours shown in Figure |5j results in an esti- 
mated filament width of 0.15 — 0.27 pc. Because this 
widths a based on the column density threshold, they 
might be larger than the corresponding FWHM o f a fit- 
ted Gaussian profile. Arzoumanian et al. (|201 1|) calcu- 

: ilam( 



tersby et al. (2011) prov ided a large scale tem perature 
and column density map, Beuther et al.| ( |2012[ ) present a 
more detailed study of IRDC 18454 adjacent to the mini- 
starburst W43. In contrast to IRDC 18454, G48 provides 
a temperature map of an externally undisturbed fila- 
ment with the better resolution of the convolution to the 
SPIRE 350 /im beam (as for G030.49+0.96, IRDC 18454 
data were convolved to the SPIRE 500 /im beam) . The 
low temperatures in the cloud centers indicate that they 
a re very young s t ructu res . 

|Wilcock et al.| ( |2012| ) used a different approach. They 
compared radiative transfer models of starless purely ex- 
ternally heated cores with the observed embedded cores 
in IRDCs in the Spitzer and Herschel bands. They found 
core temperature profiles with minimum temperatures of 
8— 11 K which increase by ~ 10— 17 K towards the edge of 
the core. While our temperature mapping only observe 
an average line-of-sight temperature of the optically thin 
material, the RT model can retrieve temperature pro- 
files along the line of sight. Since these profiles probing 
also the optically thick regime, the minimal temperatures 
are obviously lower. On the other hand, the RT models 
require an assumption about the geometry of the cores 
and cannot reassemble the diffuse temperature structure 
of the whole IRDC. 

The temperature profile for G48 shows a drop towards 
the inner part of the detected sub-millimeter peaks. The- 
oretical considerations show that gas-phase coolants, as 
C + and CO, are efficient for low densities, while the en- 
ergy transfer between gas and dust is t he dominant cool- 
ing process for highe r densities (e.g., |Goldsmith 2001 



Q- 8 P C ; jPeretto et al.| ([2012} and |Arzoumanian et al.| [Glover fc ClarkpOl^. With the presence oi IRAC anc. 



lated the line mass for their filaments. All filaments have 
line mass below the average line mass for G48 and most 
of the filaments are in contrast to G48 below the crit- 
ical thermal line mass. Also most of th e filaments in 
the P ipe Nebula are thermally subcritical ( |Peretto et al.| 
|2012[ ) . These results might indicate a trend of larger fila- 
ment widths towards supercritcal filaments (with respect 
to turbulence). 

4.5. Temperature profiles 

Recent studies have shown that the temperature struc- 
ture of star-forming regions is non-isothermal. Dust tem- 
perature maps for low-mass star-form ing Bok-globules 
are shown in, e.g., |Stutz et al.| ( 2010| CB244). Prelim- 
i nary results for IRDCs were shown by |Peretto et al. 



MIPS 24 /im point sources one would expect to see the 
imprint of this embedded heating sources in the temper- 
ature maps. The SPIRE 350 /im beam size limits the 
resolution of our temperature maps to ~ 0.3 so that 
the embedded heating sources remain unresolved. The 
dust temperature mapping as presented in this work can 
provide a versatile tool to test the predicted coupling 
between gas and dust. Methylacetylene (CH 3 C=CH) 
can be used to determine gas temperatures in massive 
star- forming regions (cf. Miettinen et al.||2006[ ). IRAM 
pointed observations of CH 3 CCH towards the peak posi- 
tion of the sub-millimeter peak PI results in a gas tem- 
perature of 26 ± 7K. Although this gas temperature is 
7K higher than the dust temperature the difference is 
not significant due to the high uncertainty in the gas 
temperature estimate. However the trend to higher gas 
temperatures is a good motivation for further observa- 



tions of CH 3 C= 



:CH or NH 3 . Especially gas temperature 



( 12010] ) for G030.49+0.96 using SED-fitting. The min- 



mapping can provide useful information on the gas-dust 
coupling at higher densities in molecular clumps. Future 
work should study the relation between dust tempera- 
ture structures an d gas temperature maps as provided 
by ammonia (as in Pillai et al.||2QQ6 Ragan et al.||2011 ) 
or methylacetylene line observations. Such data should 
provide further observational evidence for the predicted 
strong thermal coupling between gas and dust at higher 



IRDC G048.66-0.29 
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densities. 

The northern and southern extent of the IRDC (host- 
ing C2 and C3) show a much flatter temperature pro- 
file with higher temperatures compared to the central 
part of the filament. Furthermore, the column densities 
are lower towards C2 and C3. This leads to the conclu- 
sion that the gas and dust densities are lower in these 
parts of the G48. Pure isothermal cylind rical collapse 
can hardly explain this density gradient. Har tmann fc 
2QQ7[) showed by modeling the density struc 



Burkert 

tures lor the Orion A cloud that, even in the isothermal 
case, a simple initial arrangement of a slightly rotating 
gas sheet can result in filamentary structures with com- 
plex density structures. 

As in the studies mentioned above, it is a clear find- 
ing of this work that the dust temperatures are non- 
isothermally distributed in an IRDC. What impact could 
the modeling of thermal profiles have? Thermal line 
broadening is significantly smaller compared to the ob- 
served linewidths. Therefore, a change of local temper- 
atures should not have a notable influence on the tur- 
bulence structure inside the core. On the other hand 
the thermal profile towards the core can play an im- 
portant role in the support against fragmentation, since 
the Jeans mass is strongly dependent on temperature as 
Mj ~ T 3 / 2 p -1 / 2 . Theoretical studies of clustered star 
formation utilizing polytropic equations of states P ~ p 1 
and temperature dependent density profiles T ~ p 7_1 in- 
dicate that the number of fragments in a cloud is strongly 
dependent on the polytropic index 7 ( |Li et al" 2003). 
For the non-isothermal case (7 / 1) the polytropic in- 
dex reflects the balance between cooling and heating. 
Deriving dust temperatures maps as in this work is the 
first step towards a better understanding of spatial dis- 
tributed cooling and heating inside the core. A more 



detail ed discussion of t hermal properties can be f ound in 



e.g. 



Larson 



(2005) and Commercon et al. (2011) 



5. SUMMARY 



The following gives a brief summary of the results pre- 
sented in this paper: 

• G48 is a quite isolated, mono- filamentary IRDC. 
The GMC W51 observed in the same field is twice 
as distant as G48. 

• The dark filament observed by Spitzer does not ap- 
pear as a large coherent extinction silhouette in the 
PACS bands. Only the region of highest extinction 
in the central part is seen in extinction. The fila- 
ment is clearly seen in emission in the SPIRE bands 
beyond 250 /im. This is in good agreement with the 
sub- millimeter emission maps from SCUBA and 
ATLASGAL. 

• We present reliable temperature and column den- 
sity maps based on pixel-by-pixel SED fitting of 
Herschel, ATLASGAL and SCUBA legacy maps. 
Regions of low emissions outside of the galactic 
plane were utilized to estimate the global back- 
ground level for the Herschel maps. The location 
of the temperature minima are well aligned with 
the sub- millimeter and column density peaks. The 
temperatures drop to 17.5 K in the central part of 



G48. The heating by embedded MIPS 24 fim point 
sources is not detected, due to beam dilution. 

• We derived a total mass of 390 M for G48 by ap- 
plying a threshold of 1.1 • 10 22 cm -2 to the SED 
fitted column density maps. 

• "Green and fuzzy" features reveal the central parts 
of the IRDC as the sites of most active ongoing star 
formation. 

• We detected 13 point sources in all three PACS 
bands. Only one source lacks a MIPS 24 /am coun- 
terpart. 

• SED fitting to these sources result in core masses 
of 0.6 — 14.3 M , while the bolometric luminos- 
ity ranges between 15 and 630 L . The two 
most massive cores (Source 1 and 2) have mass 
of 14.3 and 11.9 M and a bolometric luminosity 
of 630 and 475 L , respectively. Comparing their 
light-to-mass ratio to evolutionary models, both 
sources are likely to evolve into high-mass stars 
with M > 8 M . Applying the definition from low- 
mass star formation, four cores are good candidates 
for "Class 0" objects. 

• The average line mass of the G48 filament is 
174M pc _1 . Compared to theoretical predictions 
the whole filament is collapsing, even under the as- 
sumption of internal turbulent support. 

• The point source separation in G48 is significantly 
smaller than in other IRDCs, but in good agree- 
ment with the predictions by the isothermal cylin- 
drical collapse model. 

• According to our virial analysis the molecular 
clumps found in G48 are gravitationally unbound 
or just pressure confined. 

• Our observati ons are consistent w ith previous ob- 
servations by Ormel et al. (2005), which deduced 
increasing turbulent motion towards the inner 
parts of the clumps. 

• We detected several molecular species, such as 
HCO+, HCN, HNC, N 2 H + , 13 CO, C 18 (and 
CH 3 CCH). The mapping observation showed 
several molecular clumps associated to the sub- 
millimeter peaks. 

• We found CO depletion towards the molecular 
clumps CI + C2 with factor of fo ~ 3.5. 

• The low abundance of CO towards the sub- 
millimeter peaks PI and P2 permit relative high 
abundances of N 2 H + . 
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